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A generalised computer programme has been used to calculate the hourly total 
cooling load for a cold storage for a given orientation The heat capacity of the 
wall is also included in terms of time-lag and decrement factor. 

Result have been obtained for the case of surface evaporation. It shows that 
surface evaporation renders about 5 % saving in cooling load. But a saving of 
about Rs 80,000 per year is found for a normal size cold storage Due to surface 
evaporation the temperature of roof is almost constant and lesser which reduces 
thermal stress in building material It has been also found that there is decrease 
in cooling load due to reduced storage temperature as at low temperature the rate 
of respiration is minimum. 

To solve energy problem in cold storage, the ammonia-absorption system is pro- 
posed operating on the primary energy. This reduces operational cost by 40 % and 
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also helps the cold storage owner to get rid of problem of frequent power break 
down 

A mathematical model is developed for the cooling load prediction. The model is 
tested for its applicability on a cold storage. The theoretical prediction compares 
very well with that of the actual cooling load of the cold storage. 
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Chapter 1 


Introduction 


1.1 Description 


India being an agricultural country produces large amount of perishable agri- 
cultural commodities of different varieties. Some 20 to 30 % of these products 
are wasted in the areas where they are grown, because of nonavailability of 
controlled environment [1]. 

Of many perishable products the preservation of potatoes is the major applica- 
tion covering 80 % of the cold storages in India [2]. The use of cold storages has 
been further diversified to store other commodities like fish and dairy products 
(5 %), vegetables, flowers, spices and other items (15 %). 

Cold storages have been proved to be boon to potato cultivators. Their growth 
in India has been quite rapid being functionally related as 

F{x) = 3145 -H 9256 4x - 66.646x2 + 92 5552x3 (1.1) 

where x = Christan era - 1947,as, 1947 has been taken as base in obtaining 
the above correlation. 

With increasing demand of cold storages for various locations, rapid growth 
in this industry is taking place. In terms of capacity of cold storages, Uttar 
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Pradesh tops the list with West Bengal at the second place. 

Cold storages maintain a low-temperature and suitable humidity to prevent 
the spoilage of perishable products, make them available in off season and also 
in other regions where they are not harvested. The growers also get due share 
of the profit by not selling their seasonal goods at throw away prices. 

The rapid developments that are taking place in the cold storage industry 
has led to greater emphasis on the design aspect. Unlike the arbitrary design 
of earlier years, new vistas has been opened up in the cold storage designs. 
These days it is most common practice to construct cold storages and refrig- 
erated warehouses with larger volumes, capable of storing a larger amount of 
commodities at a lesser cost. Improvements are being made in this field to 
minimise the functional energy requirements and the cost wherever possible. 

Generally, a number of alternatives can be thought of for designing a new 
system to bring of improvement the existing one. In cold storage the structural 
cooling load which is dependent on environment is quite substantial apart from 
commodity load. As we can not reduce commodity load which is fixed one, so 
our main aim is to reduce structure cooling load. 

In this regard, energy saving can be achieved by decreasing the solar radiation 
falling on walls and roofs of the cold storage by sinking the structure partly or 
fully below the ground level. Another advantage of sinking the structure partly 
or wholly under ground is that the underground temperature of the soil is lower 
than the ambient temperature and moreover the underground temperature 
fluctuation with time is very small as compared of ambient temperature. 

The insulation requirement for underground cold storages is considerably re- 
duced than that for cold storages 

above ground level, as a result of lesser heat transmission through underground 
walls. 


The concept of an underground cold storage however has some disadvantages; 
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like extra digging cost and higher underground constructional cost and unfore- 
seen problems due to seepage of water from the underground. 

The shape and orientation of the structure also plays a vital role in energy 
saving. Since the west and east facing walls are exposed to solar radiation 
for the longer period. As such the area of these walls should preferably be 
smaller than those of north and south facing walls. Other orientation can 
be resorted to if results in decreased structural cooling load and lesser costs 
The best shape of the cold storage is that which renders minimum surface 
area for a given volume, rendering reduced cooling load and lesser insulation 
requirement. The shape of the cold storages is influenced by other practical 
problems like the constructional cost, maintenance and handling cost etc. 

Lastly, energy saving can also be achieved by selecting proper refrigerating 
machinery and its regular maintenance in order to have a high coefficient of 
performance. Hence, the optimum cold storage design is envisaged as having 
both financial feasibility and that of functional energy which renders the overall 
cost for the cold storage to be minimum. 


1.2 Literature Review 


The work available in the field of design of refrigeration and air-conditioning 
systems, about surface evaporation/evaporative cooling and related things per- 
taining to cold storages and refrigerated warehouses is summarized in the fol- 
lowing paragraphs: 

1.2.1 Structure Heat Gain 

The analysis of cooling load involves parameters such as the ambient air tem- 
perature, direct and diffused solar radiations, air velocity, characteristic of 
enclosing walls and the orientation of the building. 
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Threlkeld [3] has extensively described the process of heat transfer due to 
temperature difference between the external and internal environments. The 
periodic heat transfer through the walls has been explicitly dealt with, using 
the concept of sol-air temperature His analysis essentially deals with heat 
transfer through homogeneous walls. 

Mackey and Wright [4] have described in detail the periodic heat transfer 
through composite walls or roofs by reducing a composite wall problem into 
an equivalent homogeneous wall. The following three methods account for the 
periodicity of external conditions: 

1 Threlkeld ’s classical approach [3] 

2 Transfer function method [5j 
3. Finite difference method [6] 

Threlkeld’s classical approach involves the concept of sol-air temperature and 
is used in the present work. 

In the transfer function method the various components of space heat gain 
are added together to get an instantaneous total rate of space heat gain. It 
is then converted into an instantaneous space cooling load through the use of 
weighting factors called coefficients of ‘room transfer function’. The transfer 
function is nothing more than a set of coefficients that relate to an output 
function at a given time to the value of one or more driving functions at the 
time of previous times. 

Kadambi and Hutchinson [6] have described an approximate technique to de- 
termine one dimensional transient heat transfer through walls and roofs, in 
the form of the Finite Difference Method. The basic simplicity of approximate 
method contrasted with analytical techniques is asserted in their work. 

A detailed procedure of cooling load calculations, from the point of view of 
practical designs, is outlined in ASHRAE Handbook of Fundamentals [7] . The 
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methodology and equations for hour-by-hour load calculations are presented 
in the handbook. 


1.2.2 Economic Design and Optimization 

With the limited availability of energy and its consequent increasing cost, it 
has been necessary to use energy optimally. Extensive work has been done 
in minimizing the functional energy and cost of cold stores and refrigerated 
warehouses which are above the ground level or partially underground. 

McClure[8] has described a method for the optimum design of building systems 
to reduce energy requirements. 

Bonar [9] has analysed the different parameters and factors that affect the 
economics of a refrigeration system for a warehouse wdth reference to the min- 
imization of total costs. Thermal insulation is found to be the foremost factor 
that affects the functional energy requirement and the energy^ cost as well. 

Hemmings [2] has provided a break-up of the total cost of cold storage as 
follows: 

building work (30%), steelwork and cladding (18%), refrigeration plant and 
electric bill (21%) and insulation (31%). Thus it can be seen that insulation 
itself accounts for one- third of the cost. This gives a statistical picture as to 
the method of controlling the total overall cost of the whole system. 

The interested article of Spiel vogel [10] demonstrate that the use of more in- 
sulation than required can increase the energy consumption and cost. 

Gupta [11] has done extensive work on automated optimum design of refriger- 
ated warehouses and airconditioned buildings. He has formulated the design 
problem as a nonlinear mathematical programming problem and has used mul- 
tidimensional optimization technique to solve it. His work presents a proba- 
bilistic optimum design procedure by considering the randomness of the input 
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parameters. He has used approximate partial derivative method to evaluate 
the cooling and heating load. 

Heinge [12] has made a comparative study of the single and multi-storeyed 
cold storage, and the dependence of insulation thickness on multi-storeyed 
constructions Single storey cold storages are suitable for places where the cost 
of land is quite cheap, whereas multi-storeyed cold storage has been found to 
be more economical provided mechanical handling is used. 

The work of Claesson and Holnquist [13] envisages the storage of cold and 
frozen food in underground mined rock caverns. Their works give an idea 
of the energy saving from underground storage. They have suggested that 
such storages can be preferred where disused mines are readily available for 
underground construction. 

Prasad et al [14], in an analytic study, have determined the expression for op- 
timum insulation thickness and sink of an underground cold storage in terms 
of various non-dimensional parameters. He has determined a quadratic rela- 
tionship between the approximate optimum sink and insulation thickness. 

1.2.3 Evaporative Cooling/Surface Evaporation 

The increasing cost of electric energy and frequent power break down have 
necessiated an economic method of cooling. The evaporative/surface cooling 
is one of such method. 

Leonardo da Vinci built a water powered evaporative cooler for the bed room 
of his patron’s wife for first time. It gave reasonable comfort conditions inside 
the room. 

The first mechanical direct evaporative cooler was developed in about 1932 
recently Jain [15] did extensive work in the field of roof surface evaporation. 
He used gunny bags to keep the roof wet at low cost . He also developed an 
electronic sensor for automatic starting of the pump for sewetting of gunny 
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bags. The same is used in the present study for cold storage. 

The characteristics of dry tropical climate are given in Jannot Yves [16]. The 
needs of economic equipment for air cooling are outlined from climatic data 
and local financial possibilities. For this climatic zone, an acceptable thermal 
comfort can be maintained all the year round by direct evaporative cooling 
of outside air. Experimental results obtained by testing a direct evaporative 
cooler built locally at a compatible cost with local incomes are given. 


1.3 Present Study 


As these days there is more concern for energy conservation, development of 
energy efficient system and utilisation of thermal energy directly is the main 
emphasis in present work to deal with energy efficient system for a long term 
solution of energy problem in cold storage. 

For the accurate evaluation of the cooling load of a cold storage a generalised 
computer programme has been developed by considering daily outdoor tem- 
perature variation based on actual as well as approximate value. The same is 
compared with actual cooling load. 

The surface evaporation on roof of a cold storage renders a saving in the 
cooling load has been found out to the tune of 5 %. The effect of lower surface 
temperature has been also studied. At the same time it will reduce the thermal 
stress in building materials. 

The cooling load is also calculated by decreasing the storage temperature. 
Decreasing the storage temperature , cooling load gets reduced significantly. 

Cold storage faces frequent power cut and suffers from frequent breakdowns 
causing reduced profit and loss of stored commodity. To solve this problem, the 
replacement of the existing ammonia vapour-compression system by ammonia- 
water vapour-absorption system is recommended. The energy and economic 
analysis reveal that by using ammonia-water vapour-absorption system, the 
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operating cost is reduced by 40 %. The investment for the new system is 
recovered in 2.5 years when 15 % compound interest is considered. 

Economy analysis based on energy cost is carried out for Lithium Bromide- 
Water Absorption system for upper side and Ammonia Vapour-Compression 
system for lower side. Unfortunately proposition flouts the very purpose of 
rendering electric power to bare minimum. 



Chapter 2 


Evaporative Cooling 


Evaporative air cooling occurs in nature near waterfalls, flowing streams over 
lakes and oceans, under summer showers and even upon wet skin. The classic 
example of evaporative cooling which we read in elementary school science 
books is the cooling effect felt, when a moistened hand is waved in the air. The 
evaporative process simply removes sensible heat (i e., cooling by decreasing 
the surface temperature) and replaces it with latent heat (i.e., increases the 
moisture content of air)[21,22]. 


2.1 Primitive Evaporative Cooling 


Evaporative cooling was probably one of the first mechanical cooling methods 
used by men. Egyptian painting from 2500 B.C. shows slaves fanning porous 
clay jars to get chilled water. A wall painting in Herculaneum of 70 A.D. 
shows a leather w'ater bottle used for cooling water by similar means. Both the 
American Indians of the South West and the ancient Perscian cooled their 
tents by damp felt or grass mats. Leonardo Da Vinci built a water powered 
evaporative cooler for the bed room of his patron’s wife. 

In the ancient period in India, evaporative cooling was even used to make ice. 
For this shallow beds dug in the earth were filled with 300 mm. straw, upon 
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which shallow earthen pans were placed. For still frosty nights, even though 
air temperature fell no lower than 6°C , ice would formed sometimes 40 mm 
thick due to adequate refrigeration produced by evaporation and radiation in 
the night sky. The people here, still use the ancient cocos tatti wdth greater 
effectiveness. Doors and window's were replaced by tatties covered with dried 
khuss-khuss grass. These were kept wet traditionally by coolies, now' replaced 
by pumps with catch basin. 


2.2 Modern Evaporative Cooling 


In modern days, evaporative coolers are extensively used in textile mills, indus- 
trial plants and homes after 1900 A.D. Textile mills were natural applicants 
as they need both cooling for comfort and high humidity to keep the fibres 
workable. By 1953, evaporative cooling had become a $ 30 million per year 
industry, providing people with the means to cope with the summer climate. 
Before automobile air-conditioning was available, small evaporative coolers 
were used in the South-West to cool cars • Even today, in India people use 
khuss-khuss tatties on the roof of the car for evaporative cooling. 

The first mechanical direct evaporative cooler w'as developed in 1932 when 
house holders assisted their wet-cloth hung wundows with electric fans. By 
1933, thousands of home made coolers had appeared in Arizona. At first, 
these were burlap covered wooden frames wet by dripping water and mounted 
at outside wundows. In 1934, burlaps were replaced by excelsior pads. 

Today evaporative coolers are commonly used in residential and commercial 
buildings in tropical region with limited commercial and industrial applica- 
tion else where. There are also application of evaporative cooling in large 
commercial HVAC systems. 
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2.3 Theory 


Evaporation is described as an adiabatic process, meaning that the total 
amount of heat in the thermal system remains constant. As water evaporates, 
the sensible heat content of the system falls, while the latent heat content in- 
creases by an equal amount. In other words, dry-bulb temperature falls, but 
moisture content rises. The limit of temperature reduction is up to the air 
wet-bulb temperature at the beginning of the process. Surface evaporation 
can cool up to the wet-bulb temperature. The process stops when the relative 
humidity approaches 100 %. A simple measure of the potential for evaporative 
cooling at any given air condition is the wet-bulb depression, or the difference 
between the dry-bulb and wet-bulb temperature It provides an upper limit of 
the achievable temperature change by direct evaporation. 

The evaporative cooling results from a material phase change, i.e., in particular 
evaporation of water. The amount of heat absorbed by a material as it changes 
from liquid to vapour is called the latent heat of vapourisation. For water, the 
latent heat of vapourisation is given approximately by: 

(2500 - 2.3 X water temperature in °C ) kJ/kg (2-1) 

The initial temperature of the water being evaporated has some effect on the 
net heat absorbed. However, for most processes, it can be assumed that each 
kilogram of water that evaporated, absorbs 2442.5 kJ of heat, corresponding 
to water temperature of 25 "C. 

A perfect evaporative cooling process is illustrated in Figure 2.1 for typical May 
condition at Kanpur. The outdoor condition ( 41.5 ‘’C dry-bulb temperature 
and 22.7 % R.H. ) is represented by point A on the psychrometric chart. Air 
drawn from building interior is at 32 °C and 40 % R.H., represented by C. 
Outside air is sent through the evaporative system and leaves it at saturation 
condition, B. This idealized process occurs along a constant enthalpy and stops 
at 100 % R.H. The humid air at B is then mixed with internal air at C. The 
final indoor air condition will lie some where along the line from B to C. Its 
exact location depends on the relative masses of the two air streams. A mixture 
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DBT — 32 41.5 


Figure 2.1: Direct evaporative cooling 

of one half evaporatively cooled air and one half interior air would provide air 
at 28 °C and 63 % R.H 


2.3.1 Actual Evaporative Cooling Process 


Actual evaporative cooling systems have inherent inefficiencies and limitation 
that prevent them from operating as ideal process. The idealized processes 
may be considered as the absolute maximum attainable performance. The 
real system will always provide less cooling i.e., higher temperature. 


In an ideal evaporative system, all the air sent through the system is com- 
pletely saturated with moisture and leaves at 100 % R.H. In actual situation 
all water is not evaporated and so air is not fully saturated. The capability of 
an evaporative cooler to cool and humidify the air is measured by saturation 


efficiency rjse as: 



( 2 . 2 ) 
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Commercially produced systems provide saturation efficiencies of about 80 %, 
with some time as low as 50 % and other as high as 90 %. 


2.4 Evaporative Cooling equipment 


Evaporative cooling equipment may use either a direct or indirect cooling 
process. In the direct process the w^ater is evaporated directly into the air 
stream that flows into the conditioned space. In an indirect cooler the evap- 
oration process is separated from the air to be delivered to the conditioned 
space. Evaporative coolers may also be either single -stage or multi-stage. 
Multi-stage evaporative systems use two evaporation processes, wdth the first 
supplying precooled air to a second stage. 

2.4.1 Direct Evaporative Cooling 

In direct evaporative cooler water is supplied through a float valve to a small 
reservoir and then flows down through fibrous pads. A fan draws large vol- 
umes of outdoor air through pads, where it is cooled by evaporation, and then 
supplied to the building. This cool and more humid air absorbs sensible heat 
from the building. 

The temperature of air delivered by an evaporative cooler may be estimated 
by the following equation. 

'^supply ~ T DBout B depression ^ ^se) {2'^) 

or, 

Tsupply TDBout - {TDB out TWBout) X Vse (2.4) 

Using the typical saturation efficiency, i.e., 80 % this equation can be rewTitten 


as 


Tsuppiy = 0.2TDBout - O.STWBout 


(2.5) 
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This equation is useful for judging the applicability of direct evaporative coolers 
under different conditions. When conditions permit their use, evaporative 
coolers are much more economical method of cooling than conventional vapor- 
compression air-conditioners. The power consumption of fan, which ranges 
from 250 W to 750 W in most residential units. A comparable air-conditioner 
might consume six to eight times as much electric power. Water consumption 
in an evaporative cooler is quite significant. A typical evaporative cooler with 
a saturation efficiency of 80 % requires about 9 liters of w'ater per hour of air 
flow for each 6 °C of wet-bulb depression. 

An evaporatively cooled building operates as an open system where dry out- 
door air is continuously drawn into cooler, conditioned, and supplied to the 
building, providing a complete air change every 2 to 4 minutes. 

Direct evaporative coolers humidify the air supplied to a building, rendering 
the relative humidity of indoors to be always higher than that of outdoors. The 
successful application of direct evaporative coolers depends on the existence of 
outdoor humidity levels well below human comfort conditions. But even in arid 
zones, they are considered as second-class cooling system, because they cannot 
humidify consistently as done by conventional air-conditioner. Consequently, 
evaporative systems should be installed with conventional air conditioning sys- 
tems An airconditioning system which incorporates evaporating cooling and 
mechanical airconditioning system is called hybrid air-conditioner 


2.4.2 Indirect Evaporative Cooling 

Indirect evaporating coolers attempt to make use of evaporative cooling process 
wuthout increasing the amount of moisture in the supplied air to the building. 
Indirect coolers use a heat exchanger to separate the direct evaporative process 
from the air to be delivered to the building. A direct evaporative process cools 
air that flow^s across one side of the heat exchanger, removing heat, and is then 
exhausted to the outdoors. The air to be supplied to the building flows across 
the other side of the heat exchanger and is cooled without receiving moisture 
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directly. 


The cooled water from the cooling tower is circulated through the heat ex- 
changer by a small pump. A fan drawls the air through the heat exchanger and 
supplies it to the building. The system is quite effective in dry climates, often 
reducing indoor temperature by 10 to 20 °C on a hot and dry climate. A sim- 
ple indirect cooling process is shown on the psychrometric chart in Figure 2.2. 
Outdoor air at 41.5®C dry-bulb temperature and 22.7 % R.H. enters the cooler 
at point A and cooled without the addition or removal of moisture along a hor- 
izontal line to point B, 25°C dry-bulb temperature and 57 % R.H. (the actual 
output conditions depend on the saturation efficiency of the direct evaporative 
process and the effectiveness of the heat exchanger). 
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Figure 2.2; Indirect evaporative cooling. 


Two-stage evaporative cooling system offers significant improvements in per- 
formance over indirect evaporative cooling systems. In a two-stage evaporative 
system, the cooled air from an indirect cooler is passed through a direct evap- 
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orative cooler for additional cooling. 


2.5 Roof Surface Evaporation 


Roof surface evaporation [15], a simple and passive process of cooling buildings 
is effective, economically viable and practicable technology to solve all the 
global problems of cooling of buildings. This process not only eliminate major 
flow of heat to the building through its roof, but it also extracts heat from the 
building as shown in fig 2 3. Evidently, the surface evaporation brings the roof 
temperature close to the wet-bulb temperature. Hence it becomes a heat sink 
and so heat gain by the building through the walls, windows etc. is transferred 
to the roof, rendering reasonable comfort temperature inside the building. This 
type of cooling is quite suitable for workshops, industrial buildings, theaters 
etc 

Figure 2.4 shows a schematic diagram of the roof surface evaporation. It 
comprises a gunny bag layer over the roof, water sprinkler having a pipe con- 
nections, a pump to maintain adequate pressure in the sprinkler and a water 
tank. The evaporative cooling needs small quantity of water upto 9 liters per 
square meter of roof surface area per day under severe hot-dry weather condi- 
tions. Water has to be sprayed uniformly over a thin absorbing material lining 
providing wick action on the entire roof surface. It must absorb sufficient wa- 
ter and evaporate it quickly. The heat and water absorbing material lining 
is required to be maintained moist day and night during summer periods for 
continuous and quick evaporation. 

The past manual spraying on empty cement bags for low cost housing schemes 
has been replaced by highly refined system to meet the requirement of sophis- 
ticated modern buildings. This has led to development of automatic spraying 
device and coirmating of the surface as heat absorptive material lining Such 
arrangement is reliable and renders better performance. It is the most eco- 
nomic cooling system. Since evaporation takes place at all temperatures and 
based on well known relationship of latent heat of water versus its tempera- 




Figure 2.3: Flow of heat towards roof. 


ture, it has been found after calculation that 3800 liters of water can produce 
nearly 900 tons of refrigeration effect if it is properly used in this process as 
per its basic requirements. Quality of water do not make much difference for 
this cooling process except replacement of metallic sensing probes in automatic 
spraying system. 

Recent development of hand-operated device and its successful implementation 
at various places, has made possible for easy adoption of this process by masses 
even in remote areas where both electricity and tap w'ater may not be available. 

The precautions, which have to be taken care of, particularly for water-proofing 
aspects at the design stage by changing the topography of the roofs more- 
sloppy pitched or cylindrical dome to avoid stagnation of water on roofs. 
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Figure 2.4: A typical arrangement of surface evaporation. 

The use of conventional heat insulating materials ,e.g., lime concrete, mud 
phaska, thermocole having so many practical problem can be straight way 
eliminated once this process is made integral part of buildings. 

In the present work Roof Surface Evaporation is proposed for cold storage. 




Chapter 3 


Estimation of Cooling Load 


3.1 Different Heat Sources : 


The heat load in cold storages arises from the following sources: 

1 Energy transfer from the surroundings to the system These are: 

(a) Heat transmission through barriers such as walls, doors, ceiling, 
floor, etc. being caused by the temperature difference existing on 
the two sides of the barrier. 

(b) The solar heat gain: 

Absorbed by walls or roofs exposed to radiation from the sun and 
transferred to the inside space. 

(c) Heat and moisture introduced to the conditioned space through 
infiltration and ventilation air. 

2. Heat generation within the confined space. They are: 

(a) Heat load from the product stored inside - sensible and latent 

(b) Heat load from fans, lights and other electrical appliances in use. 
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3.2 Need for Hourly Cooling Load Calcula- 
tions 


It IS well known fact that the ambient air temperature varies during the day 
and its ^a^iation also depends upon the day of a month. Accordingly, the 
heat load varies over a period of 24 hours This is mainly due to the variation 
in the solar intensity falling on the earth surface. The load calculation is 
further complicated by the fact that a wall has thermal capacity, due to which 
a certain amount of heat passing through it is stored and is transmitted to 
the inside at some time latter. Therefore, calculation based on instantaneous 
heat transmission through structure without considering the thermal capacity 
of the wall, is erroneous. Moreover, the traditional methods of evaluating the 
cooling load involves various assumptions like the load from each component 
is constant in a day and maximum value of the cooling load is equal to the 
individual maximum, which is not correct. Since the occurrence of maximum 
load at particular time is not the occurrences of the other components. It is 
proper to adopt an hourly cooling load calculation to facilitate an economic 
selection of a refrigeration system. 


3.3 Heat Transfer Through Structure: 

3.3.1 Building Survey 

The building survey will include the flowing particulars; 

1 Location of the building ,i.e., longitude and latitude 
2. Orientation 

3 Dimension of the building structure, such as length, breadth , height and 
thickness of each layer of the building material. 

4. Composition of building material and their physical properties. 
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3.3.2 Hourly Outside Temperature Variation 

Figure 3.1 shows the hourly variation in outside temperature for Kanpur. The 
metrological data reveals that the minimum temperature occurs just before 
sunshine (may be one or two hours) while the maximum temperature occurs 
three to four hours after the solar noon [3]. Since, the temperature time 
relation is available for only a few places, an hourly variation in temperature 
is predicted based on maximum and minimum temperatures. It has been taken 
in the following form: 

T = L + M cos{lbt - N) (3 1) 

where, T=temperature at any time in °C. 

L, M and N are constants, calculated by applying the boundary conditions- 

rfT' 

— = 0,forT = Tmax or Train 
at 

It. is assumed that occurs one hour before the sunshine and Tmax occurs 
12 hours thereafter. 

The comparison of the temperature variation by equation 3.1 with actual vari- 
ation is also shown in figure 3.1. The data has been taken from the average 
temiierature of three years i.e. 1993, 1994 and 1995 [18] 

3.3.3 Solar Radiation 

Solar radiation [7] forms the significant part of cooling load in the cold stor- 
ages The total radiation (/«), reaching on terrestrial surface is the sum of the 
direct solar radiation {hr), the diffused sky radiation {hf) and the solar radi- 
ation reflected from the surrounding surface (h)- The intensity of the direct 
component is the product of the direct normal radiation (hn) and the cosine 
of the incidence angle {6) between the incoming solar rays and a line normal 
to the surface as shown in Figure 3.2. Thus, 
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Figure 3.1: Hourly variation of actual and theoretical temperature. 




3 3 Heat Transfer Through Structure: 



Horizontal Surface 

Figure 3.2: Direct solar radiation on horizontal surface. 

It — Idn COS 9 + Id} + Ir (3-2) 

In the pr('sent analysis the reflected radiation (1^) is neglected. 

I, In and l,i} are given as: 

exp{B/sinP) 

Idf = C.Idn.Fss (3.4) 

value of /%,,=(). 5 for vertical surfaces, 1.0 for horizontal surfaces and (1+cos 4>)/2 
for any other inclined surfaces. 

The value of A, B and C are tabulated in [7] and are given in Appendix ‘B’. 

3.3.4 Determination of Incident Angles 


The sun’s position in sky is most conveniently expressed in terms of solar 
altitude 3, above the horizontal and solar azimuth angle 'y, measured from 
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south I'hoso angle f3 and 7 depends on the latitude of the place I, solar 
(lech nation r/, and hour angle h P and ') is related with I, d and h by the 
following expression; 


sviP = cosi .cosd .cosh + stnl.sind 


(3 5) 


and, 


COS7 = 


sinP.sinl — sind 


cosP.cosl 

'J'li(‘ houi angle is calculated as follows; 


LST = LCT + (Equation of time) 


(3.6) 


(3.7) 


Wheie, LST is local solar time and LCT is local civil time which are based 
on the standaid rneiidian of the country. And due to the irregularities of the 
eai til's rotation, obliquity of earth’s orbit and other factors, a solar day is not 
exactly ecjual to 24 hours. To counter that effect equation of time is added in 
IX ’T, which is given by 

LCT = [1ST - (82 5 - longitude)] /I'b (3 8) 

Where IS'I is Indian Standard Time which is based on 82.5 meridian (standard 
longitude foi India) passing through Naini near Allahabad 

'I'he hour angle is zero at solar noon and increases by 15® every hour. Hence 
hour angle, 

h — (12 - LST) X 15, before solar noon (3.9) 

h = (12 + LST) X 15, after solar noon (3.10) 

If a suiface is tilted by an angle (f) to the vertical position then the incidence 
angle 0 is given by 

(•osB = cosP.cosoi.cos4> -h sinp.sincj) (3-11) 


If till' surface is vertical, then 


cos6 = cosp.cosoi 


(3.12) 



3 3 Heat Transfer Through Struc ture : 

If the surface is horizontal, (^= 7 r/ 2 , then 


cos9 = sin^ 


(3.13) 


Aft('r knowing the value of h, the values of /? and 7 are calculated The wall 
solar azimuth angle, a [7] for various walls are- 


Ok = 

Oh' = 


K/2 - 7 I 

k /2 + 7l 

, before noon 
, after noon 

^ east facing wall 

(3.14) 

l7r/2-7| 

i7r/2-k7| 

, before noon 
, after noon 

^ west facing wall 

(3.15) 

K-tI , 

k -t- 7i , 

before noon \ 
after noon j 

north facing wall 

(3.16) 

as = 

7 , south facing w'all 

(3.17) 


3.3.5 Sol-air Temperature 

The calculation of space cooling load as a result of heat gain through exterior 
roof and walls involves the concept of sol-air temperature [3]. A heat balance 
at a sunlit surface gives the heat flux to the surface as; 

q = arit + ho{To - T^,o) - eAi? (3 18) 

(= 1 , for black body 

A/f=0.0G3 kW/m^, for horizontal surface and 0, for vertical surface [7] 

The above equation can also be written as 

q = hoiTsoi - T^,o) (3-19) 

where, 


Tso/ = To Oir h/ho - c i^Rjho 
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3.3.6 Heat Transfer Through Walls and Roof Using Decre- 
ment Factor and Time Lag 


As building materials have finite thermal capacity, the instantaneous heat 
transfer through walls and roofs by neglecting thermal capacity of the building 
materials is not correct The structural detail of walls and ceiling is given in 
figure 3 3 and 3.4 respectively. 



GLASS WOOL 


K1 =0 0007 kW/mC 
K2 = 0 000039 kW/mC 

Figure 3.3: Structural detail of wall 
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Figure 3.4: Structural detail of ceiling. 
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Thermal capacity of wall Q/,, is expressed as' 

Cth = me = peV = pcAAt (3.20) 

Due to thermal capacity of a structure the following two effects have been 
observed [17] 

1 . There is a time lag between the heat transfer at the outside surface and 
the inside surface. 

2. There is decrement in heat transfer due to the absorption of heat by the 
wall and subsequent transfer of a part of this heat back to the outside 
air when temperature of ambient is lower. 

Considering these two facts, actual heat transfer at any time ’t’ is given by: 

Qt = jziU.A.iZa - T,) + (3.21) 

; = 1 

Q, = AU{Ti - Z) (3.22) 

where. Qb is the basement load and Z is the underground temperature being 
taken as 25°C. 

Hence. 

Qs=Qt + Qb (3-23) 

The value of 01 and A with respect to wall thickness are given in Figure 3.5 
and 3.6 respectively. 

In the present study equation 3.21 is used to calculate the heat transfer through 
the structures. 

3.3.7 Over-All Heat Transfer Coefficient for Walls and 
Roof 

To get overall heat transfer coefficient U , the value of ho and h, are taken from 
[6], as given below: 


0.2 
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Figure 3.6: Variation of decrement factor with wail thickness. 
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For exterior wall 

/?o = 0 00391 (0.399 + t^,) (3.24) 

For roof 

/i, = 0.00141(1.454 + 1;^) (3.25) 

For the inner surfaces of the enclosed space, the convective heat transfer coeffi- 
cient IS found to be the function of the inside air velocity and the temperature 
difference existing between the room air and inner surface of the wall. In terms 
of interior air velocity v, and temperature difference ’A Tj, we get 

For interior walls 

/i. = Cl ( AT0° + 0.00391 (0.399 + u,) 

Foi ceiling 

/i, = C2(A7;)° 0.00141(1.454 

where, f'l and are constants and its values are 

Cl = 1 25 X 10"^ C2 = 1.45 X 10'^ 


Now, the value of U for walls and roof taking into account of the time-lag is 
calculated as follows. 


g = U[{Zoi-Z) + XTait-^i)]=^UD 

(3 28) 

where, 

(3.29) 

D = {Tsoi — Ti) -1- X{To(t-^l) — Tern) 

For the interior wall: 


q = [ci (AT, )° + 0.00391(0.399 + v^)]AT 

(3.30) 

q = ci(ATl)^'^^ + h-ujATr 

(3.31) 


(3 26) 
(3.27) 


w’here 


= 0.00391(0.399 + ^;^) 
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Comparing the equation 3.28 and 3.30 we get 




But 


Putting 


U = 


f + E^ + f 

t%o A.J tlx 


1 _ ^ 

H— \- — = R = Constant 

K-t rif) 


(3 32) 

(3.33) 

(3.34) 


for a given wall and for constant outside wind velocity(uu,) equation 3.32 be- 
comes 


Cl (AT,)' + AT. = - 


D 


D 


1 25 


+ R 


D = AT. + i?Ca(AT.)' + RKAT^ 


(3 35) 
(3 36) 


The Newton-Raphson iterative procedure has been used to get AT. Accord- 
ingly we have 

F = AT. + Fci(AT.)' 25 + Rh^AT^ - D (3.37) 

Thus, after getting value of AT., h, is calculated and hence values of U'‘ for 
different walls and roof are calculated 


3.4 Infiltration Load 


There is alw^ays leakage of outdoor air into a building through cracks and 
openings caused by pressure difference across the boundary surfaces. This 
leakage is called infiltration load The exchange of air may cause both types 
of load namely latent heat load and sensible heat load. The load due to 
infiltration is given by the following expression: 


Qtnjil 


^ Toom^ach{,Ro,a Ri,a) 

r>„,.(24 X 3600) 


(3 38) 
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3.5 Lighting and Other Electrical Application 


In general the instantaneous rate of heat gain from electric lighting is calculated 
from 


Q light = Total Light Wattage x Cooling Load Factor (CLF) (3.39) 

CLF is a function of time of use, type of arrangement, room furnishing etc in 
present study, a value of 0.88 is taken from[7] 

Other electrical appliance used normally is fan and the instantaneous load 
from the fan is equal to the total rating((5/an) 

Therefore, 

Q equip — Qlight Qf an (3.40) 

3.6 Sensible and Respiration Load of Stored 
Products 

In cold storage generally potato is stored. Hence, 

• _TN X Resp 
24x 3600" 

During loading period, there is sensible heat load which is equal to 

Qsen = TN X 1000 X Cpp{To - T) 

3.7 Total Load 


(3.41) 

(3.42) 


Total load is sum of all the load, i.e. Solar load, infiltration load, equipment 
load, respiration load and sensible load. 

Q total ~ Q solar ~t“ Qinftl “b Q equip + Qr "b Q sen (3 43) 



3.8 Solar Load With Surface Evaporation 


32 


3.8 Solar Load With Surface Evaporation 

When surface evaporation is applied to the roof of the building, the tempera- 
ture of the roof surface approaches to wet-bulb temperature- Hence the solar 
load on the roof with the surface evaporation is given by 

Qroof = AnUR{T^,t-T.) ( 3 - 44 ) 

Using the procedure outlined in this chapter, the hourly cooling load and hence 
total yearly cooling load are calculated with and without surface evaporation 
for "Arnar Cold Storage” Chaubepur, Kanpur and the same is compared with 
the actual load. 


Chapter 4 


Ammonia- Absorption System a 
Substitute of Vapour- 
Compression System for Cold 
Storage 

4.1 Introduction 

Cold storage plays significant role in preventing the spoilage of the huge quan- 
tity of perishable products by providing an appropriate environment. Thus, it 
serves the nation in a big way and also helps grower to sell their produce at 
premium in place of selling the same at throw away prices. 

Unfortunately, these days cold storages are facing a serious problem of short- 
age of power, especially in rural areas. The existing system depends fully 
upon external power supply on which they have no control. So this calls for 
development of a system which may operate at disposal of the owner of a cold 
storage. In the present case ammonia- absorption system is studied to replace 
ammonia-compression system. Analyses are carried out on the basis of data 
collected from local cold storages. It is found that operational cost of the 
ammonia-absorption system is about 60 % to that of existing system. The 
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detailed result are given in Chapter 5. The compressor of the existing system 
is replaced by a generator, an absorber, an analyser, a rectifier and a pump as 
evident from a typical ammonia-absorption system shown in figure 4.1. By this 
arrangement major components of the compression refrigeration system and 
their orientation and installation remain unaltered. This aspect convinces cold 
storage owner for accepting the change. Despite high initial cost compared to 
the existing system, the additional cost is compensated within five to six years 
due to saving in the running cost. 



EXISTING COMPONENTS 


Figure 4.1: A typical ammonia-absorption system for a cold storage. 


4.2 Why Cold Storage Should Shift to Vapour- 
Absorption System ? 

As cold storage is necessary for perishable products,the existing erratic power 
supply leads to spoilage of the commodities. Hence, a self dependent system is 
the required solution. Moreover, the existing power scenario reveals the dismal 
position in continuous or uninterrupted power supply to the cold storage. This 
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leads to procurement of a big Diesel-Generating set at an enormous investment 
(about 3 to 4 lacs). This additional standby unit raises the cost of preservation 
load on diesel fuel requirement. 

A cold storage generally stores 6000 to 7500 tonnes of potatoes having 100 
kW of power consumption. If there are 10 cold storages, the total powder 
requirement becomes 1000 kW, i.e., 1.0 MW. To supply 1.0 MW of power at 
the cold storage, the installed capacity of the power house should be 1.0 x 
1.25/0.6 = 2.0 MW as 25 % of pow'er is lost in transmission and plant load 
factor of power plant is about 60 %. Hence, to get 1.0 MW at the place of 
use the installed capacity of 2.0 MW requires an investment of Rs. 9 crores. 
It implies, that the Government has to invest Rs. 90 lacs to meet the energy 
need of each of cold storages. The net result is the increased cost of potato 
storage. To reduce the storage cost, the theft of power or tempering of meter 
reading by some owners cannot be ruled out. 

Hence, adoption of ammonia-absorption refrigeration system would not only 
solve the power problem of cold storage owners, but also help to improve power 
situation for commercial, industrial and domestic activities Moreover it would 
be beneficial to both the owmer of cold storage and the Government. 


4.3 Existing and Modified System 


Figure 4.2 exhibits the existing system and the suggested ammonia- absorption 
system with common components It is evident that only compressor is getting 
replaced by a generator, an absorber, a pump and a heat exchanger. The rest 
system remains unaltered. Hence, there is no need of any change inside the 
cold storage. So the owners may get convinced for the change in their existing 
system. 





TO BE REMOVED 
WHEN SYSTEN CHANGED 


Figure 4.2- A comparision between existing and change needed in the existing 
system. 









Chapter 5 


Result and Discussion 


This chapter deals with the calculation of cooling loads on the basis of the the- 
oretical and actual temperature variation having effect of heat capacity and 
time-lag factor of structures For this a generalised computer programme is 
made based on mathematical formulation discussed in chapter 3. Flow dia- 
gram is shown in figure 5.1. To have energy conservation surface evaporation is 
applied on the roof of cold storage which reduces the roof temperature at one 
side and nullify the effect of radiation on other side. It is also discussed that de- 
crease in inside temperature of cold storage, there is reduction in cooling load 
due to decreased respiration load of potatoes. To solve energy problem am- 
monia vapour-compression is replaced by ammonia-absorption system which 
reduces operating cost substantially. 


5.1 Graphical Representation of Results 

5.1.1 Theoretical and Actual Temperature Variation 


The actual hourly temperature variation as well as the suggested sinusoidal 
hourly variation in temperature based on maximum and minimum temperature 
variation of the day have been shown in figure 3.1. The actual temperature is 
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Figure 5.1: Flow chart. 
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generally seen to be slightly higher than the theoretical value from 8 A.M. to 
2:30 P.M., just before peak hour. Thereafter, the actual temperature is lower 
than that of the theoretical value. The areas under both the curve reveal a 
small difference (3 to 8 %). 

5.1.2 Sol-air Temperature Variation For Different Walls 
and Roof 

The ambient temperature and sol-air temperatures over the different walls and 
the roof is presented in figure 5.2. The typical result presented is for the month 
of May. From this plot it is seen that the sol-air temperature for a horizontal 
surface (roof) is the highest and is maximum at noon, corresponding to the 
maximum incident solar radiation From the plot it can be inferred that, the 
heat transfer through the east and west facing walls is more prominent than 
that of south facing walls. Thus the area of eastern and western walls should 
be less than that of northern and southern walls. 

Figure 5.3 shows the hourly variation in wet- bulb and temperature of different 
walls and roof for a particular day. The wet-bulb temperature is much lower 
than that of temperatures of different walls and remains almost constant. This 
vindicates that the evaporation in nature is close to an adiabatic saturation 
process. Figure 5.4 shows temperature variation of roof and wet-bulb tempera- 
ture. If surface evaporation is applied on the roof, the roof surface temperature 
tends to w'et-bulb temperature. Thus, we see that there is substantial reduc- 
tion in roof surface temperature and hence reduced cooling load by applying 
surface evaporation Figure 5.5 exhibits variation in dry-bulb temperature and 
wet-bulb temperature for a typical month of May. 
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Time(hrs) 


Figure 5.2: Variation of ambient temperature and sol-air temperature for walls 
and roof for the month of may. 

5.1.3 Variation of Wet-Bulb Temperature and Relative 
Humidity 

Figure 5.6 shows variation of relative humidity with time for the month of May. 
The relative humidity is minimum at about 3 P.M. when dry-bulb temperature 
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Figure 5.3: Variation of ambient temperature and sol-air temperature for walls 
and roof along with wet-bulb temperature for the month of may 


is maximum. The nature of the curve is the same for each month. This helps 
to keep the wet-bulb temperature almost constant (variation range 3 tom 5 °C 
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) as shown in figure 5 7. The wet-bulb temperature variation is maximum for 
the month of November, December and January as during that time maximum 
minimum temperature and relative humidity variation is maximum. Where as 
wet-bulb temperature is almost constant for the months of July and August 
when maximum minimum temperature variation is minimum. 

5.1.4 Solar Load With Surface Evaporation 

Figure 5.8 shows hourly variation of solar load with and without surface evap- 
oration where as figure 5.9 shows monthly variation of solar load with and 
without surface evaporation. From the figure it is clear that hourly variation 
of solar load is almost constant with surface evaporation. This is supported 
from figure 5.7 where the hourly variation in wet-bulb temperature is found 
to be very small as compared to the ambient dry-bulb temperature.However, 
when there is no surface evaporation, it is maximum at around 8 A.M. due 
time lag factor of about 19 hours. Solar load is maximum for the month of 
June in both with and without surface evaporation. 


5.2 Tabulated Results 

Table 5.1 shows average solar load per day with and without surface evapora- 
tion, respiration load per day for each month and total load with and without 
surface evaporation per month. Generally, in cold storage, loading starts after 
15^'' February and it is filled to half of its capacity during that month. So, 
in February total load is calculated for 14 days and for half of the capacity. 
Generally, cold storage is full by the end of March but in some cases it may 
go also in first week of April. In these calculations it is assumed that, cold 
storage is full by the end of March. During February and March sensible load 
of potato is also included in respiration load. Inside temperature is maintained 
at 2 °C. After March there is no sensible heat load. Hence respiration load is 


constants. 



R H(%) 
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Figure 5.6: Hourly variation in relative humidity for may. 


Table 5.2 shows total cooling load with and without surface evaporation, per- 
centage saving in cooling load and saving in terms of rupees. The cost of 
equipment for surface evaporation is about Rs. 75,000 which will be recovered 
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Figure 5.7; Hourly variation of wet-bulb temperature for different months. 


within one year. 

Table 5.3 shows the reduction in cooling load as storage temperature decreases. 
This is for the month of May. This is due to reduced heat of respiration. Heat 
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Figure 5.8' Hourly variation of solar load in may. 



Figure 5.9: Monthly variation of solar load. 
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of respiration is taken from [20]. The same is given at temperature 4.5 °C and 
0 °C. Linear interpolation is used for heat of respiration at other temperatures. 

Table 5.4 shows actual cooling load supplied by Amar Cold Storage, Chaubepur, 
Kanpur and cooling load calculated using computer programme. The data sup- 
plied by cold storage is meter reading, i.e., energy consumed by compressor. 
The same is compared with calculated load by converting in terms of cooling 
load by multiplying with COP. 

Actual cooling load of February is much less than the calculated load. This 
is due to fact that in calculated load it is assumed that half of loading is 
completed m February but the same is not actually true. In calculation it is 
assumed that loading is over by the end of March. But actually loading has 
gone to the month of April. That is the reason for more actual cooling load 
than the calculated load in the month of April. For the other month cooling 
load IS comparable 


Table 5.1: Solar, respiration and total load. 


Month 

Solar Load 

Respiration 

Load 

Total Load 

Without 

Evap. 

With Evap. 

Without 

Evap. 

With Evap. 


kWh/day 

kWh/day 

kWh/day 

kWh/month 

kWh/month 

Feb. 

592.22 

449.59 

2222.22 

42216.56 
for 15 days 

40077.71 
for 15 days 

March 

774.84 

591.69 

4444.44 

156578.63 

151084.1 

April 

972.42 

744.8 

3413.33 

131572.64 

124744.02 

May 

1091.21 


3413.33 

135136.33 

127945.34 

.Tunc 

1126.29 

IH^HI 

3413.33 

136188.81 

129359.55 


1011.9 

[SEBH 

3413.33 

132576.92 

127509.51 



817.2 

3413.33 

131714.17 

126911.11 

Sept. 

949.99 


3413.13 

130899.64 

126274.48 

IliiM 

IBIBl 

EQBHHI 

2666.67 

106292.48 

101496.24 

Nov. 


587.82 

1706.67 

72940.27 

68834.48 
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Table 5.2‘ Cooling load with and without surface evaporation and saving. 


Total cooling load without surface evaporation 

1176296.4 kWh 

Total cooling load with surface evaporation 

1124241.5 kWh 

Saving in cooling load 

52054.9 kWh 

Percentage of saving 

4.42 % 

Saving in terms of Rs. 

Rs. 78082.32 


Table 5.3: Reduction in cooling load due to reduced storage temperature. 


Storage 

Heat of 

Total Load 

Temperature 

Respiration 

Without Evap. 

With Evap. 

° C 

kj/day-ton 

kWh/month 

kWh/month 

4 5 

2720 

179334.36 

168167.00 

4.0 

2556.7 

167027.7 

159859.89 

30 

2230 

150409.1 

143240.41 

20 

1903 

135136.33 

127945.3 

1.0 

1579 7 

117177.52 

110006.98 

0.0 

1250 

100559 09 

93389.69 


Table 5.4: Actual and calculated cooling load. 


Month 

Calculated Load(kWh/month) 

Actual Load(kWh/ 

month) 

cop 3.0 

cop 3.25 

cop 3.5 

Feb. 

42216.56 

16920 

18330 

19740 

March 

156578.63 

131577 

142541.25 

153506.5 

April 

131572.64 

164325 

178018.75 

191712.5 

May 

135136.33 

131316 

141259 

153202 

June 

136188.81 

121101 

131192.75 

141284.5 

July 

132576.92 

119151 

129680.25 

139009.5 

Aug. 

131714.17 

116835 

126571.25 

136307.5 

Sept 

130899 64 

122301 

132492.75 

142684.5 

Oct. 

106292.48 

97389 

105504.75 

113620.5 

Nov. 

72940.27 

65055 

70476.25 

75897.5 

■■i 

1176276.39 

1085970 

1176467.8 

1266965 
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Table 5.5 shows the saving in operational cost due to adoption of ammonia- 
absorption system in place of vapour-compression system. This calculation is 
made from the energy data supplied by local cold storages which is given in 
Appendix D. COP of the existing vapour-compression system is taken as 3.0 
where as COP of proposed ammonia -absorption system is taken as 0.6 based 
on primary energy. Primary energy used for calculation is coal. Calorific value 
of coal is taken as 16,000 kJ/kg, and its cost is Rs. 1500 per tonne. The cost 
of power is taken as Rs 3.75/kWh as supplied by cold storage owner. 


Table 5.5: Saving in operational cost in ammonia-absorption system against 
vapour-compression system. 


Kanaoj Cold Storage, Kanauj 

Storing Capacity 

Electric Bill 

Cost of Coal 

Saving 

Percentage Saving 

7000 tonnes 

Rs. 1447500/year 
Rs. 900000/year 
Rs. 547000/year 
37.8 % 

Amar Cold Storage Chaubepur, Kanpur 

Storing Capacity 

Electric Bill 

Cost of Coal 

Saving 

Percentage Saving 

6400 tonnes 

Rs. 1357470/year 
Rs. 825000/year 
Rs. 532470/year 
39.25 % 

Raj Cold Storage, Kanpur 

Storing Capacity 

Electric Bill 

Cost of Coal 

Saving 

Percentage Saving 

11000 tonnes * 

Rs. 3053780/year 
Rs 1875000/year 
Rs. 1178780/year 
38.6 % 


* This cold storage has an ice plant also. 


Chapter 6 


Conclusions and Suggestions 


6.1 Conclusions 

1 The total cooling load has been calculated on the basis of hourly tem- 
perature variation based on Tmax and T^tn for a given day. The same 
compared with actual cooling load based on actual hourly temperature 
variation. 

2 The overall heat transfer coefficient has been calculated accurately by 
considering the dependence of inside convective heat transfer on the in- 
side surface temperature and air motion. Newton- Raphson method has 
been used to calculate the same. 

3. Surface evaporation is applied to the roof of the cold storage. It renders 
reduced cooling load to the tune of 4.5 %. This percentage is found to 
be low due to dominance of the respiration load of potatoes. 

4 By decreasing the storage temperature, the cooling gets reduced due to 
reduced heat of respiration of potatoes at lower temperature. 

5. To solve the power problem in cold storage, the replacement of the exist- 
ing ammonia-compression system by ammonia-water vapour-absorption 
system is recommended. This would enable the use of primary energy 
for refrigeration. 
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0 By using aininonia-absorption system, the operating cost is reduced bv 
40 %. 

7 . Tiie recoininended system will operate at the disposal of the owner of 
cold storage'. Hence, better quality of the stored products. 


6.2 Scope for Future Work 

1 . The surface evaporation should be applied to the roof and actual reduc- 
l.ion in cooling load should be measured. 

2. liicononiic analysis can be done including the cost of surface evaporation 
and absorption system. 

3. Aminonia-absorptioii system is to be designed and the same is installed 
in cold storage in order to demonstrate its salient features. 


Ct 
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Appendix A 


Temperature and Relative 
Humidity of Outside Air 


Table A.l: Average temperature and relative humidity of outside air at Kan- 
pur, January to June. 


mmm 

January 

February 

March 

April 

May 

Ju 

ne 

IHHH' 

T 



^ _ 

T 

4> 

lI 




T 

<f> 


" C 


UH 

wmm 

® C 




H91' 


■ISBi 

% 


10 28 

86 2 


81 45 

18 65 

B»i|| 


52 8 

28 05 

52 25 

■jb 

68 5 

5 30 

9 5 

88 9 

11 4 

85 3 

16 65 



59 4 


58 8 

WL 

74 0 

8 30 

1 1 35 

87 0 

14 7 

79 75 

21 1 

59 0 

27 8 

41 4 

32 8 

41 4 


63 1 

11 30 

18.75 

58 5 

23.0 

46 0 

29 1 

35 3 

35 15 

23 8 

39 2 

28 4 


48 5 

14 30 

22 05 

45.15 

25 7 

35 3 

33.65 

26 3 

37 25 


41 5 

22 7 

E WM 

40 25 

17 30 

19 9 

55 

23 85 

41 9 

30 65 


35 3 

23 75 

39 6 

27 25 

E 

44 3 

20 30 

13 65 

77.1 1 

17 15 

67 3 

23 95 


28 55 

38.55 

32 6 

41 3 

E ^9 

56 5 

23 30 

11 5 

83 15 

14 7 

76 5 


63 0 

25 35 

45.95 

29 78 

48 0 


64 0 
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Table A. 2; Average temperature and relative humidity of outside air at Kan- 
pur, July to December. 




August 

September 

October 

November 

December 

mmi 

r 


T 


T 

<P 

T 

4> 

T 


T 

<i> 

warn 

^ C 


O Q 

% 

^ C 

% 

® C 

% 

^ C 

% 

^ C 

% 

2 30 

27 5 

90 0 

26 85 

93 55 

25 25 

89 45 

20 3 

85 7 

14 3 

84 65 

10 7 


5 30 

26 95 

91 95 

26 55 

93 5 

24 3 

91 8 

18 95 

88 9 

12 7 

89 3 

9 2 


8 30 

29 0 

84 55 

28 2 

89 55 

27 51 

72 35 

24 3 

72 9 

17 75 

79 85 

13 1 


1 1 30 

32 0 

72 4 

30 8 

80 15 

31 25 

65 75 

31 8 

44 15 

27 2 

41 35 

22.4 


14 30 

33 2 

68 25 

31 2 

76 3 

32 5 

60 3 

33 2 

40 05 

29 3 

32 65 

24 65 


17 30 

32 45 

71 85 

30 0 

81 25 

30 5 

68 68 

28 55 

57 15 

23 3 

59 15 

20 3 


20 30 

29 6 

81 85 

28 15 

89 

27 25 

82 6 

23 5 

74 05 

18 2 

76 4 

13 3 


23 30 

28 4 

86 1 

27 35 

91 5 

26 05 

86 7 

21 05 

81 75 

15 85 

82 5 

12 05 



Table A 3: Maximum and minimum temperature of outside air at Kanpur. 


Month 

Maximum Temperature 

Minimum Temperature 


" C 

° C 

January 

23.1 

8.0 

February 

26.5 

10.05 

March 

32.65 

15.4 

April 

38.15 

19.9 

May 

41.9 

25.18 

June 

41.35 

27.6 

July 

34.55 

26.25 

August 

32.62 

26.04 

September 

33.6 

23.85 

October 

34.36 

18.23 

November 

30.15 

11.75 

December 

26.05 

8.1 












Appendix B 


Values for Solar Intensity 
Calculation 


Month 

Equation 
of Time 
minute 

Declination 

Angle 

degree 

to 

B 

C 

Jan 

-11.2 

-20 

1230 

0.142 

0.058 

Fob 

-13.9 

-10.8 

1214 

0.144 

0.600 

March 

-7.5 

0.0 

1185 

0.156 

0.071 

April 

1.1 

11.6 

1135 

0.180 

0.097 

May 

3.3 

20.0 

1103 

0.196 

0.121 

June 

-1.4 

23.45 

1088 

0.205 

0.134 

July 

-6.2 

20.6 

1085 

0.207 

0.136 

August 

-2.4 

12.3 

1107 

0.201 

0.122 

Sept. 

7.5 

00 

1151 

0.177 

0.092 

Oct. 

15.4 

-10.6 

1192 

0.160 

0.073 

Nov. 

13.8 

-19.8 

1220 

0.149 

0.063 

Doc. 

1.6 

-23.45 

1233 

0.142 

0.057 


A = Apparent solar irradiation at air mass = 0, W/m^ 
B=- Atmospheric extinction coefficient 
C= Diffuse radiation factor 



Appendix C 

Cold Storage Particulars for 
Cooling Load Calculation 


Cool Room Size 
Aioa of East Facing Wall 
Area of West Facing Wall 
Area of North Facing Wall 
Area of South Facing Wall 
Area of Roof 
Area of Base 
Thickness of Wall 
Thickness of Roof 
Capacity of Cold Storage 
Number of Fans 
Rating of Fans 


37.8X 31.0 X 16.8 
31.0 X 16.8 w?. 
31.0 X 16.8 m^. 
37.8 X 16.8 
31.0 X 16.8 m2. 
37.8 X 31.0 m2. 
37.8 X 31.0 m2. 
0.56 m. 

0.60 m. 

6400 tonnes. 

84. 

40 watts. 



Appendix D 


Energy Data of Cold storages 


Table D.l. Power taken by compressors of Kanauj Cold Storage, Kanauj of 
capacity 7000 tonnes 


Month 

Compressors Power 


kWh 

Feb. 

22000 

March 

47000 

April 

47000 

May 

52000 

June 

57000 

July 

47000 

August 

52000 

Sept. 

32000 

Oct. 

17000 

Nov. 

13000 
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Table D.2: Power taken by compressors of Amar Cold Storage, Chaubepur, 
Kanpur of capacity 6400 tonnes 


Month 

Compressors Power 


kWh 

Feb. 

5640 

March 

43859 

April 

54775 

May 

43772 

June 

49367 

July 

39717 

August 

38948 

Sept. 

40767 

Oct. 

32463 

Nov. 

21685 


Table D.3: Power taken by compressors of Raj Cold Storage, Kanpur of ca- 
I)acity 11000 tonnes 


Month 

Compressors Power 


kWh 

Jan. 

15325 

Feb. 

15690 

March 

81943 

April 

134437 

May 

127266 

June 

116033 

July 

75467 

August 

86637 

Sept. 

80850 

Oct. 

45810 

Nov. 

21453 

Dec. 

13430 





Appendix E 


Storage Condition and 
Properties of Potato 


Storage Temperature 
Inside Relative Humidity 
Freezing Point 

Sp. Heat Above Freezing Point 
Sp. Heat Below Freezing Point 
Heat of Respiration 


2 ° C. 

85 %. 

-1.7 ° C. 

3.6 kJ/kg ° C. 

1.97 kJ/kg " C. 

1250 kJ/day-ton (at 0 " C) 
2720 kJ/day-ton (at 4.5 ° C) 



Appendix F 


Properties of Moist Air 


221 287 

Pi ~ cx):)[{7 21379+{1 152x10-^-4 787xa0-®7'<i6)(Td6-4S3 15)2}{647 31/r<i6-l}] 


Pv = Pvh 


(Pa ^u'h){'^db '^vjb) 

1940-1 UT^b 


BclativeHumidity{(p) = ^ 



